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摘要 : 近 50 年 来 ,气候 变化 加 剧 、 全 球 变 暧 明 显 ,陆地 表面 趋向 干旱 ,与 此 同时 ,森林 病 虫 灾害 损失 不 断 增 加 ,两 者 关系 密切 。 以 
江西 省 为 研究 区 域 ,采集 1961 一 2010 年 气温 .湿度 .日照 等 气候 要 素 及 1992 一 2010 年 森林 病 虫 灾害 发 生 面积 程度 数据 ,获知 
研究 区 过 去 50 年 年 均 气 温 .冬季 均 温 都 在 以 0.16C/10a、0.27%C/10a 速率 在 上 升 ,而 年 均 相对 湿度 则 以 -0.45%/10a 的 速率 在 
下 降 , 显 示 其 气候 变 暖 .干旱 化 趋势 显著 。 同 时 ,其 森林 病 虫 灾害 面积 则 以 58125hm*”/10a 速率 扩 增 。 研 究 以 相关 分 析 、 主 成 分 
分 析 法 对 所 选 24 个 气候 要 素 降 维 筛选 ,以 逐步 回归 法 构建 模型 ,以 小 波 分 析 研 究 区 内 时 空 变化 等 分 析 两 者 关系 。 得 到 16( 与 
病害 发 生 ) 或 17( 与 病虫害 与 虫害 发 生 ) 个 显著 相关 气候 要 素 ,最 大 正 相 关 要 素 为 夏季 均 温 滑动 均值 等 ,而 最 大 负 相 关 要 素 三 者 
均 为 温 湿 系数 滑动 均值 ( 温 湿 系 数 = 年 均 相 对 湿度 / 年 均 温 ) ;获得 4 或 5 个 主 成 分 ,而 代表 温度 或 温 湿度 联合 变量 特征 主 成 分 
贡献 率 最 高 (病虫害 .虫害 :41.43% ,病害 :42% ) ;建立 3 个 以 森林 病 虫 灾害 为 纵 轴 ` 温 湿 系 数 滑动 均值 作 横 轴 且 具 预测 能 力 的 
最 优 回归 模型 (了 ws =3.582 x 105-7.750 x 105X、J = 一 6.375x105X+2.95x105 与 Ya =-1.375x105X+6.321x105) ,其 线性 拟 


合 度 分 别 为 77.9% .79.1% 与 56.7% ,平均 预测 准确 率 分 别 为 66.2% .68.6% .47.9% , 而 研究 区 温 湿 系 数 滑动 均值 过 去 50 年 在 显 
著 下 降 , 并 在 1993 年 后 转向 负 距 平 。 这 说 明 气 候 变 暖 与 环境 干旱 化 对 区 内 森林 病 虫 灾害 发 生 具 驱动 作用 ,促使 其 发 生 加 重 , 尤 
其 在 进入 20 世纪 90 年 代 后 表现 更 烈 。 得 知 温 湿 系数 滑动 均值 具 29a 周期 变化 , 若 以 1993 年 为 起 始 ,不 考虑 非 气候 要 素 作 用 ， 
则 其 病虫害 加 重 态势 将 持续 到 2022 年 。 在 空间 上 ,过 去 50 年 闵 南 温 湿 系 数 滑动 均值 最 低 , 预 判 其 森林 病 虫 灾 害 将 较 其 他 区 域 
发 生 更 重 ,而 闵 东 , 北 温 湿 系数 滑动 均值 倾向 率 最 高 , 预 判 其 森林 病 虫 灾害 发 生 的 变动 性 可 能 更 大 。 以 上 结果 说 明 ,气候 变 暖 与 
环境 干旱 加 重 江西 省 森林 病 虫 灾害 ,六 南 应 为 重点 防治 工作 区 域 ,要 加 强 闵 东北 的 监测 预警 工作 。 
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Abstract: Over the past 50 years, the climate has changed rapidly with global warming and land surface drying, which has 
been accompanied by increased forest loss and damage resulting from pests and diseases. Climate data (e. g., air 


temperature ，relative humidity, and sunshine data) from 1961 to 2010 showed that the climate in the Jiangxi Province was 
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warming significantly, with the annual mean temperature increasing by 0.16°C per 10 years, and the winter mean 
temperature rising by 0.27°C per 10 years. The Jiangxi Province climate was also drying throughout this period (annual 
mean relative humidity decreased —0.45% per 10 years ). In addition, forest pests and disease occurrence from 1992 to 
2010 showed that in Jiangxi, the area affected by of forest diseases and pests increased significantly, with 58,125 hm” per 
10 years. Pearson correlation and principal component analyses showed that 16 (for forest diseases) or 17 (for forest pests ) 
climate elements were significantly related to the forest loss. From these individual elements ，the most positively correlated 
was a 9-year sliding average of summer mean temperature, and the most negatively correlated component was a 9-year 
sliding average of hydrothermal coefficient (annual mean temperature / annual mean relative humidity ) . Amongst the four or 
five principal components, the variables temperature and temperature-humidity contributed most to explaining forest area 
loss (41.43% and 42.0%, respectively ). In stepwise regression analyses, three optimal regression models ( Total : 了 = 
3.582 x 106 -7.750 x 10°X, forest pest: 了 = -6.375 x 10°X + 2.95 x 10°, forest diseases: Y=—1.375 x 10°X + 6.321 x 
10’”) were analyzed to describe the forest loss (Y) by the 9-year sliding average of temperature humidity coefficient (X). 
The three models showed a linear fit of 77.9%, 79.1%, and 56.7% and a prediction accuracy of 66.2%, 68.6%, and 
47.9%, respectively. A declining trend in the sliding average of temperature humidity coefficient was observed over the past 
50 years, for which the anomaly transferred from positive to negative in 1993. This indicates that climate warming and 
droughts could have aggravated the forest loss and damage over the past 50 years, especially after 1990s. A wavelet analysis 
showed a 29-year periodicity in the temperature humidity coefficient. If this anomaly started in 1993, the forest loss trend 
could potentially be relieved by the end of 2022. In Jiangxi, the temperature humidity coefficient significantly increases from 
the south to the north, suggesting that forest diseases and pest disasters should be more extreme in Gan Nan than in other 
areas. In addition, since the change rate of climate warming or drought was higher in Gan Dong and Gan Bei, a high 
variability of forest diseases and pest disasters can be expected there in the future. Overall, our results suggest that climate 
warming and environmental drought aggravates forest diseases and pest disasters in Jiangxi. They furthermore emphasize that 
Gan Nan could be a key area in preventing and controlling the effects of forest diseases and pest disasters, whereas the 


monitoring efforts in Gan Dong and Gan Bei should be increased. 


Key Words: long-time; global warming; droughts; forest loss generated by pests and diseases; driving effect; regionality 


20 世纪 以 来 ,全 球 地 表 升温 加 快 ,过 去 50 年 升温 速率 是 过 去 100 年 的 2 倍 '" ,20 世纪 90 年 代 是 近 百 年 
来 最 暧 的 10 年 "1, 且 过 去 30 年 的 每 个 10a 地 表 平 均 温度 都 高 于 1850 年 以 来 的 任何 一 个 10a 等 1。 伴随 气 
候 变 化 加 剧 全球 变 暖 明 显 、 及 陆地 表面 趋向 干旱 化 5 等 ,各 类 灾害 频繁 发 生 ( 如 森林 灾害 '“* ) .风险 上 升 . 损 
失 加 重 * 中。 和 森林 生物 灾害 是 破坏 森林 生态 系统 ,降低 “地 球 绿 肺 ”功能 的 重要 因素 ,在 气候 变 暧 驱动 下 ,其 
破坏 力 在 增强 ' ,但 仍 待 提供 更 多 证 据 。 据 范 霍 夫 定律 ,气候 变 暖 可 促使 生物 发 育 速率 加 快 . 生 殖 力 增强 . 存 
活 率 尤其 是 越冬 存活 率 提高 ,从 而 提高 次 年 发 生 基数 水 平 。 气 候 变 上 暧 也 可 促使 有 害 生 物 适 生 区 变 宽 ,发 生 始 
见 期 \ 迁 飞 期 高 峰 期 提前 ,年 生育 期 延长 .发 生 代数 增加 ,爆发 周期 缩短 ,种 间 关 系 动荡 复杂 等 ,从 而 使 生物 灾 
害 复杂 多 变 程度 加 重 '1 ,损失 因而 旦 递增 趋势 '* 中 I。 研究 多 年 气候 变化 与 生物 灾害 关系 ,有 助 于 了 解 气 候 
变化 背景 下 生物 灾害 发 生 规律 ,依据 关键 要 素 预 测 未 来 灾害 动态 ,可 为 有 害 生物 预警 与 防 控 提 供 重 要 参考 。 
然而 ,森林 地 理 分 布 及 其 生态 系统 服务 多 呈 地 带 性 区 域 性 特征 性 5 , 需 选 定 一 个 区 域 ,研究 其 内 气候 变化 对 
森林 生 而 物 灾害 的 影响 ,获得 具 一 定 代表 性 的 研究 案例 。 

江西 省 地 处 亚热带 季风 气候 区 ,森林 生态 系统 发 育 良好 , 履 盖 率 高 。 区 内 气候 变化 大 ,气象 灾害 多 且 发 生 
广泛 ,如 夏秋 干旱 .极端 低温 等 发 生 频 繁 !9 。1958 一 2009 年 区 内 年 均 气温 以 0.16C/10a 速率 升 高 "1, 同 时， 
森林 病 虫 灾害 面积 近 些 年 也 以 每 年 4440.173hm 在 增加 。 一 些 研 究 认为 年 均 气 温 上 升 ,尤其 是 暖冬 ,会 导致 
森林 病 虫 灾害 加 重 !5 ,如 杨 干 象 .美国 白 蛾 . 松 纵 坑 切 梢 小 过 等 森林 害虫 灾害 或 因 1 月 份 暖 温 而 加 重 5”) 。 而 
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一 些 害虫 爆发 可 能 由 某 些 气候 要 素 引 发 ,如 在 干旱 胁迫 与 气候 干 热 作 用 下 ,意大利 阿尔 卑 斯 山地 区 的 云 杉 叶 
蜂 易 爆 发 成 灾 '” , 锈 色 粒 肩 天 牛 也 可 能 因 年 份 降水 量 少 温 湿 系 数 低 等 因素 而 爆发 成 灾 '”" ,马尾 松 毛虫 发 生 
则 受降 水 "3 相对 湿度 .最 大 风速 影响 更 甚 , 且 据 此 构建 预测 模型 准确 率 可 达 75%21 。 显 然 ,各 类 病虫害 对 
气候 变化 的 响应 并 非 相 同 ,这 增加 了 二 者 关系 的 复杂 性 。 江 西 省 有 472 种 森林 生物 灾害 ,其 中 易 造 成 危害 者 
70 余 种 ,严重 危害 者 20 余 种 ,包括 松 材 线虫 病 、 萧 氏 松 茎 象 . 马 尾 松 毛虫 .思茅 松 毛 虫 、 黄 浓 竹 蝗 \ 松 褐 天 和牛、 
星 天 牛 . 粒 肩 天 牛 ` 云 斑 白 条 天 牛 .华人 竹 毒 峨 、 刚 人 竹 毒 峨 、 松 针 褐 斑 病 油茶 炭 痊 病 \ 杨 树 锈病 、 加 拿 大 一 枝 黄 花 
等 。 另 外 ,当地 一 些 次 生性 有 害 生 物种 类 近 些 年 的 发 生 点 、 面 也 在 变化 ,在 部 分 地 段 爆 发 成 灾 , 而 少数 偶发 性 
种 类 也 在 部 分 区 域 频 繁 发 生 并 成 灾 。 总 体 上 ,该 区 森林 病 虫 灾害 正 严重 影响 着 林业 发 展 , 如 2010 年 其 森林 生 
物 灾害 面积 57586hm? ,直接 经 济 损失 8416 万 元 ,生态 服务 价值 损失 34570 万 元 "i 。 因 此 , 摸 清 气 候 变化 背景 
下 森林 病 虫 灾害 的 变化 规律 ,无疑 将 有 利于 森林 保护 工作 的 科学 实施 ,确保 森林 生态 系统 健康 发 展 。 目 前 , 较 
多 研究 专注 于 温度 .降雨 .相对 湿度 等 单一 气候 要 素 变 化 对 森林 病 虫 灾害 的 影响 ,而 针对 气候 的 长 期 综合 作用 
的 研究 较 少 。 为 此 ,本 文通 过 采集 江西 省 过 去 多 年 的 气候 综合 变化 (如 温度 、 相 对 湿度 、 降 水 量 等 要 素 均值 、 
变 率 等 数据 ) 与 森林 病 虫 灾害 发 生 数据 ,揭示 前 者 对 后 者 发 生发 展 的 影响 及 诱 变 规律 ,为 当地 及 周边 省 份 林 
业 病 虫 灾害 防 控 工 作 提 供 参 考 。 


1 研究 方法 


1.1 研究 区 域 

以 江西 省 为 研究 区 域 , 即 24"29' 一 30?04'N ,113"34' 一 118"28 中 E 间 。 其 地 处 典型 亚热带 季风 气候 区 ,气象 
灾害 种 类 多 ,频率 高 i ,为 南方 林业 重点 区 域 之 一 ,森林 病 虫 灾 害 频繁 ,年 均 发 生 30 万 hm 左右 ,最 高 可 达 50 
万 hm。 
1.2 数据 采集 整理 

采集 研究 区 内 南昌 吉安 庐山 .景德镇 .玉山 、 郡 
阳 、 修 水 贵 溪 、 樟 树 南城. 宜春. 广 昌 、 遂 川 .赣州 和 寻 
乌 15 个 气象 站 点 (图 1)1961 一 2010 年 间 的 降水 .气温 、 
相对 湿度 日照 时 数 等 地 面 气候 资料 日 值 数据 ( 自 中 
气象 数据 中 心 http://www. escience. gov. cn/metdata/ 
page/index.html,2014.11 登录 ) , 从 全 国 、 江 西林 业 统 计 
年 鉴 或 资料 中 采集 江西 省 森林 病 虫 灾害 1992 一 2010 年 
间 发 生 面积 (包括 病虫害 病害 .虫害 面积 ) 及 1998 一 
2010 年 间 轻 中、 重度 面积 数据 。 对 于 降水 类 数据 , 先 
对 漏 测 或 降水 值 进行 检查 ,以 最 小 二 乘法 对 缺失 值 进行 
插值 处 理 后 ,修正 固态 降水 值 。 同 样 地 ,分 析 前 对 其 他 
数据 也 都 用 最 小 二 乘法 对 缺失 值 进行 插值 处 理 。 以 上 
所 选 气象 指标 , 既 考 虑 其 每 年 的 均匀 变化 (X 均值 ,图 
2) ,也 考虑 其 短期 的 突变 (了 ~- 了, ,图 2) ,同时 ,采用 
滑动 滤波 法 上 进行 平滑 处 理 ,以 确保 代表 研究 区 内 多 ! 
年 的 气候 状态 。 从 而 ,使 得 各 项 气象 指标 在 均值 滑动 图 1 江西 省 气象 资料 日 值 数据 采集 站 点 
均值 . 变 率 多 个 层面 整理 分 析 。 而 且 , 为 确保 滑动 均值 Fig.1 Weather data sampling sites of Jiangxi Province 
以 时 间 序 列 中 项 加 入 到 时 间 坐 标 ,采用 9a(9 年 ) 滑动 
平均 方式 (10a 为 1 个 年 代 , 但 滑动 均值 年 数 通 常 为 奇数 ) 来 显示 气候 要 素 均 值 . 均 方差 , 变 差 系数 的 变化 趋 
势 。 滑 动 均值 : 


。 气象 站 点 
一 省 级 行政 区 
一 市 级 行政 区 
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式 中 ,XX 表示 气候 要 素 ,t 表示 年 份 ,tt4 表示 上 年份 向 前 
滑动 4a,1-4 表示 上 年份 向 后 滑动 4a。 
1.3 气候 要 素 确 定 与 筛选 

气候 变化 是 各 类 气候 要 素 均 值 或 离 差 变化 的 综合 ， 2 
即 气 候 基 本 状态 在 变化 。 如 图 2 所 示 ,均值 向 左 偏 移 则 二 = > 


Xl 2 
状态 下 降 , 向 右 偏 移 则 上 升 ; 而 某 气 候 要 素 变 差 系数 
(或 变 率 ) 的 变化 幅度 天 小 ( 变 幅 ) “(图 2) 则 意味 着 极 均值 (与 杰 率 (六 ，_J。 ) 示 意图 


端 气候 事件 发 生 频 率 、 强 度 的 变化 '*”1。 本 文选 择 温 。 Fig.2 Sketch map of mean(X) and variability( Y,,, —Yi ) 


度 .相对 湿度 .降水 日照 等 气候 要 素 ,包括 年 平均 气温 、 
夏季 均 温 (6 一 8 月 ) 冬季 均 温 (12 月 至 次 年 2 月 ) ,温差 (夏季 均 温 -冬季 均 温 ) .年 降水 量 、 年 均 相 对 湿度 , 温 
湿 系数 (年 均 相 对 湿度 / 年 均 温 ) .年 均 日 照 时 数 均值 .9a 滑动 变 率 、9a 滑动 均值 24 个 气候 要 素 多 年 数据 。 
为 避免 共 线性 干扰 , 先 将 所 有 气候 要 素 与 森林 病 虫 灾害 发 生 面积 进行 相关 分 析 后 ,筛选 出 与 之 显著 相关 者 进 
行 主 成 分 分 析 ,之 后 选择 各 主 成 分 中 载荷 值 最 大 者 为 目 变 量 .森林 病 虫 灾 害 发 生 面积 作 因 变量 ,以 逐步 回归 法 
构建 预测 模型 ,最 终 获 得 影响 森林 病 虫 灾 害 的 重要 气候 要 素 。 
1.4 气候 要 素 时 空 特征 分 析 

在 线性 倾向 分 析 中 ,用 y, 表 示 样 本 量 为 n 的 某 一 气候 要 素 ,用 1 表示 y; 所 对 应 的 年 份 ,建立 关于 x; .ti 的 一 
元 线性 回归 方程 y=ati+6b (i=1,2,…,n;a 为 回归 系数 ,6 为 回归 常数 ) ,用 最 小 二 乘法 估计 a4 和 4 , 即 &= 


n 1 n n 

之 yi 人 冯 站 并 pn 1 n 1 n 

a a 下 - , i ,t= i; ;aX10 表示 气候 倾向 率 , 用 于 定量 描述 气候 
Se, Sie D1) i=1 i=1 
i=1 nn i=1 


序列 的 变化 趋势 .特征 ,a 为 正 数 表示 气候 要 素 y 呈 上 升 趋势 ,反之 则 旺 下 降 趋势 。 男 外 ,用 气候 异常 指数 来 
评价 气候 要 素 是 否 超出 正常 范围 | ,以 年 距 平 大 于 或 小 于 1.5 个 标准 差 表 示 异 常 偏 高 或 异常 偏 低 , 确 定 异 常 
现象 出 现 的 时 间 。 

为 把 握 研 究 区 气候 时 空 变化 特征 \ 规 律 ,对 森林 病 虫 灾害 轻 . 重 灾区 及 趋势 进行 预 判 ,以 回归 分 析 气 候 倾 
向 率 ,并 基于 Matlab .AreGIS 软件 ,用 小 波 分 析 以 上 相关 气候 要 素 的 周期 性 及 其 空间 变化 动态 。 数 据 统 计 分 
析 均 用 SPSS 17.0 软件 完成 。 


2 结果 分 析 


2.1 气候 对 森林 病 虫 灾 害 的 影响 与 趋势 
2.1.1 森林 病 虫 灾害 发 生 趋势 

由 表 1 可 知 ,研究 区 森林 病 虫 灾害 在 过 去 多 年 间 处 于 明显 上 升 态 势 , 虫 害 .病害 及 病虫害 面积 不 断 扩大 ， 
线性 趋势 显著 ( 表 1) 。 例 如 ,区 内 森林 病 虫 灾害 面积 在 1992 一 2010 年 间 平 均 以 58,125hm2Z10a 的 速度 增加 。 
然而 ,新 增 受害 面积 多 为 病 虫 灾 害 轻 度 或 中 度 发 生 区 , 而 重度 受害 面 的 增长 趋势 不 明显 。 


http://www.ecologica.cn 


2 期 张 鹏 起 ”等 :气候 变 暖 .干旱 加 重 江西 省 森林 病 虫 灾害 643 


表 1 1992/1998 一 2010 年 江西 省 森林 病 虫 灾害 发 生 趋势 与 线性 方程 
Table 1 Linear equation and trend of forest diseases and pest of Jiangxi Province during 1992/1998 to 2010 


发 生 Occurrence 年 份 (X) 线性 方程 ; 趋势 


(y,hm? 次 ) Year Linear Equation Trend 
病虫害 总 发 生 Total 1992 一 2010 7y=2.168x104X -4.315x107 0.825 站 
Diseases and Pest 轻 度 发 生 Mild level 1998 一 2010 7= 1.996x104X -3.9811x107 0.777 fe 
中 度 发 生 Moderate level 1998 一 2010 Y=2.846x103X -5.646x105 0.441 从 党 
重度 发 生 Serious level 1998 一 2010 7=2.501x103X -4.985x105 0.300 一 
病害 Diseases 总 发 生 Total 1992 一 2010 Y=3.474x103X -6.917x105 0.662 个 
轻 度 发 生 Mild level 1998 一 2010 Y=3.258x103X -6.493x10° 0.531 修志 
中 度 发 生 Moderate level 1998 一 2010 Y=—100.863X + 2.073x105 0.011 一 
重度 发 生 Serious level 1998 一 2010 Y=—238.027X -4.750x105 0.181 一 
虫害 Pest 总 发 生 Total 1992 一 2010 Y=1.820x104X -3.623x107 0.828 I 
轻 度 发 生 Mild level 1998 一 2010 7= 1.670x104X -3.331x107 0.798 全 
中 度 发 生 Moderate level 1998 一 2010 Y=1.947x103X -5.853x10° 0.502 全 
重度 发 生 Serious level 1998 一 2010 Y=2.262x103X -4.510x10° 0.284 一 
* P<0.05; * * P<0.001;-: 不 显著 ; 1 :上 升 ; | :下 降 


2.1.2 ”气候 要 素 变化 对 森林 病 虫 灾害 的 影响 

研究 区 气候 变 暧 及 干旱 化 趋势 明显 ,如 过 去 50 年 (1961 一 2010) 其 年 均 气 温 .冬季 均 温 分 别 以 0.16%C/10a 
(了 =0.0174X-34.454,R? =0.3239,P<0.001 ,7 为 年 均 温 均值 ,X 为 年 份 ) .0.27% /10a( 页 =0.0283X-30.7366 ， 
R?*=0.152,P<0.001,Y 为 冬季 均 温 均 值 ,X 为 年 份 ) 速率 在 升 高 ,而 年 均 相 对 湿度 则 以 -0.45%/10a( 了 = 
-0.0429X+ 85.199, R* =0.09646, P=0.028,Y 为 相对 湿度 均值 ,X 为 年 份 ) 的 速率 在 下 降 。 并 且 ,1992 一 2010 
数据 ( 表 2) 也 显示 年 均 温 均值 .年 均 温 滑动 均值 .夏季 均 温 滑动 均值 .夏季 均 温 均值 等 显著 上 升 ,而 相对 湿度 
均值 年 降水 量 滑动 均值 相对 湿度 滑动 均值 等 则 显著 下 降 。 其 他 如 日 照 时 数 滑 动 均值 .年 降水 量变 率 、 冬 季 
均 温 变 率 、 相 对 湿度 变 率 .温差 变 率 . 温 湿 系数 变 率 等 气候 要 素 也 在 显著 上 升 ( 表 2) , 温 湿 系 数 均值 . 温 湿 系数 
滑动 均值 则 在 显著 下 降 。 


表 2 气候 要 素 变化 及 其 与 森林 病 虫 灾害 Pearson 相关 性 
Table 2 Trends of climate elements and Pearson correlation between them and forest diseases and pest 
相关 系数 Correlation coefficient(7) 


气候 要 素 趋 病虫害 
ee a 人 i and 虫害 病害 
a Pest Disease 
夏季 均 温 滑动 均值 Summer mean temperature sliding average 0.622 从 0.838™™ 0.837 一 0.647 一 
日 照 时 数 滑 动 均值 Sunshine duration sliding average 0.515 1 0.813** 0.839** 0.612”* 
年 降水 量变 率 Annual precipitation variability 0.761 人 0.808 ** 0.835** 0.609* 
冬季 均 温 变 率 Winter mean temperature variability 0.925 作 二 0.805 ** 0.805** 0.709** 
年 均 温 滑动 均值 Annual mean temperature sliding average 0.979 Ws 0.799** 0.800** 0.700** 
相对 湿度 变 率 Relative humidity variability 0.777 | 0.742 0.769 0.547” 
温差 变 率 Temperature difference variability 0.866 但 和 0.718** 0.733** 0.576”* 
夏季 均 温 均值 Mean summer average temperature 0.393 i 0.636™* 0.614 一 0.674 一 
温 湿 系数 变 率 Hydrothermal coefficient variability 0.741 作 汪 过 0.601 0.629” 0.421 
年 均 温 均值 Annual mean temperature average 0.511 os 0.586*"* 0.568™ 0.616**™ 
温差 均值 Temperature difference average 0.050 一 0.37 0.336 0.491” 
日 照 时 数 变 率 Sunshine duration variability 0.608 人 0.569” 0.604 0.362 
相对 湿度 均值 mean relative humility 0.396 1 一 0.621 一 0.632 一 一 0.506” 
年 均 温 变 率 Annual mean temperature variability 0.162 -0.630” —0.646** 一 0.494 
温 湿 系数 均值 Hydrothermal coefficient average 0.623 上 -0.716** —0.714™* —0.654** 
年 降水 量 滑 动 均 值 Annual precipitation sliding average 0.461 小 一 0.763 ** 一 0.789 -0.572* 
相对 湿度 滑动 均值 Relative humidity sliding average 0.687 站 一 0.877 一 0.890 一 0.727 
温 湿 系数 滑动 均值 了 Hydrothermal coefficient sliding average 0.926 本 站 一 0.883 ** 一 0.890 一 0.7S3 ** 


x* P<0.05, * x* P<0.001,- 不 显著 , 1 上升, | 下降 
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不 同 气候 要 素 对 森林 病 虫 灾害 影响 各 异 ,24 个 气候 要 素 中 ,找到 17( 病虫害 与 虫害 ) 或 16( 病害 ) 个 与 森 
林 病 虫 灾害 发 生 显著 相关 要 素 ( 表 2) ,其 中 ,显著 正 相 关 要 素 包 括 夏季 均 温 滑动 均值 日照 时 数 滑 动 均值 .年 
降水 量变 率 冬季 均 温 变 率 .年 均 温 滑动 均值 .相对 湿度 变 率 夏季 均 温 均 值 . 温 湿 系数 变 率 、 年 均 温 均值 温差 
变 率 ( 与 病虫害 .虫害 ) 日照 时 数 变 率 ( 与 病虫害 ` 虫 害 ) .温差 均值 ( 与 病害 ) ; 而 温 湿 系数 滑动 均值 .相对 温 
度 滑 动 均值 .年 降水 量 滑动 均值 . 温 湿 系数 均值 .年 均 温 变 率 .相对 湿度 均值 则 为 显著 负 相 关 要 素 。 其 中 ,与 森 
林 病 虫害 .虫害 .病害 发 生 的 最 大 正 相 关 要 素 分 别 为 夏季 均 温 滑动 均值 (r=0.838) 日照 时 数 滑动 均值 (7= 
0.839) .冬季 均 温 变 率 (r= 0.709) ,而 三 者 的 最 大 负 相 关 要 素 均 为 温 湿 系 数 滑动 均值 (病虫害 :r= -0.883 , 虫 
害 :r=-0.890 ,病害 :r=-0.753) ,可 见 , 温 湿 系 数 滑动 均值 等 在 负 向 上 可 能 对 研究 区 内 和 森林 病 虫 灾害 发 生 及 趋 
势 形成 起 着 更 关键 的 作用 。 
2.1.3 影响 森林 病 虫 灾害 的 关键 气候 要 素 及 其 影响 趋势 

用 主 成 分 分 析 法 ,对 以 上 与 病虫害 .虫害 发 生 显著 相关 (17 个 , 表 2) 及 与 病害 发 生 显 著 相 关 (16 个 , 表 2) 
的 要 素 进 一 步 合 并 .筛选 ,前 者 获得 4 个 主要 成 分 ( 表 3, 即 页. 丈 、 态 . 束 ) ,其 特征 根 贡献 率 分别 为 41.43% 、 
13.9% 13.66% .7.85% , 累积 贡献 率 为 76.84%(KMO = 0.647, R=16359.97, df=136, P<0.001)。 主 成 分 了 主 
要 包含 温 湿 系数 滑动 均值 .年 均 温 滑动 均值 .年 均 温 均值 . 温 湿 系 数 均 值 夏季 均 温 滑动 均值 .夏季 均 温 均值 、 
冬季 均 温 变 率 , 年 均 温 变 率 、 降 水 滑动 均值 ,主要 体现 温度 或 温 湿 度 联合 变量 特征 ;站 ,主要 包含 相对 湿度 滑动 
均值 相对 湿度 均值 ,体现 湿度 变量 特征 ;Y, 主 要 包含 日 照 时 数 变 率 、 温 湿 系 数 变 率 、 相 对 湿度 变 率 , 体 现 气 候 
要 素 变 率 变量 特征 ;Y, 主 要 包含 降水 量变 率 、 日 照 时 数 滑动 均值 ,体现 日 照 变 量 特征 。 而 后 者 则 获得 5 个 主要 
成 分 ( 表 3, 即 了 7、、Y) ,它们 的 特征 根 贡 献 率 分 别 为 42.00%、20.81% .8.04% .6.91% .6.37% ,累积 贡 
献 率 为 84.13%(KMO =0.660, 针 =16286.514, df=120,P<0.001)。 与 上 面相 似 ,其 Y 也 主要 体现 温度 或 温 湿 


度 联 合 变 量 特征 ,Y, 主 要 体现 湿度 变量 特征 ,Y, 主 要 体现 气候 要 素 变 率 变量 特征 ,Y 主 要 体现 降水 量 大 小 年 变 
化 与 日 照 变 量 特征 ,而 7; 因 只 包含 温差 均值 , 故 代表 温差 变量 特征 。 


表 3 四 /五 主 成 分 与 气候 要 素 荷 载 矩 阵 
Table 3 Four / five principal components and matrix of climate elements 
四 主 成 分 Four Principal components 五 主 成 分 Five Principal components 
六 Y, Y; Ys 六 冯 Y; Ys 了 


气候 要 素 Climate elements 


温 湿 系 数 滑动 均值 
Hydrothermal coefficient sliding average 
年 均 温 滑动 均值 

Annual mean temperature sliding average 
年 均 温 均 值 

Annual mean temperature average 

温 湿 系 数 均值 

Hydrothermal coefficient average 

夏季 均 温 滑动 均值 

Summer mean temperature sliding average 
夏季 均 温 均值 

Mean summer average temperature 
冬季 均 温 变 率 

Winter mean temperature variability 
年 均 温 变 率 

Annual mean temperature variability 

年 降水 量 滑 动 均值 
Annual precipitation sliding average 
相对 湿度 滑动 均值 


Relative humidity sliding average 


相对 湿度 均值 


Mean relative humility 


0.980 0.101 -0.040 0.025 -0.983 0.105 -0.052 -0.006 0.059 


-0.974 -0.005 0.070 —0.064 0.978 -0.010 0.065 -0.014 -0.132 


-0.969 -0.014 0.064 -0.071 0.975 -0.029 0.050 -0.018 -0.160 


0.964 0.160 -0.026 0.053 -0.969 0.183 -0.027 0.026 0.069 


-0.951 -0.034 -0.027 0.004 0.953 -0.002 0.012 -0.005 0.210 


-0.922 -0.078 -0.021 -0.010 0.921 -0.050 0.014 -0.031 0.309 


0.843 -0.150 0.214 0.084 -0.852 -0.055 0.303 0.061 0.070 


0.638 -0.077 0.218 -0.195 


0.573 0.092 -0.134 -0.014 -0.565 -0.152 -0.433 0.056 -0.169 


0.073 0.874 -0.174 -0.155 -0.073 0.866 -0.259 -0.151 0.008 


0.067 0.856 -0.076 -0.008 -0.068 0.901 -0.125 0.023 -0.044 
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四 主 成 分 Four Principal components 五 主 成 分 Five Principal components 
气候 要 素 Climate elements 
六 Y, Ys Ys 次 y, Ys Ys Ys 
日 照 时 数 变 率 
et 0.012 0.310 0.786 0.075 

Sunshine duration variability 
泪 浊 至 数 亦 谈 

人 a a 0.191 -0.451 0.743 0.245 -0.194 -0.351 0.720 0.373 0.148 
Hydrothermal coefficient variability 
相对 湿度 变 率 


-0.118 -0.475 0. 


Relative humidity variability 


温差 变 率 


-0.054 ”一 0.328 0. 


Temperature difference variability 


年 降水 量变 率 


Annual precipitation variability 


日 照 时 数 滑动 均值 


Sunshine duration sliding average 


温差 均值 


0.054 0.079 0. 


0.005 -0.385 -0. 


661 0.281 0.091 -0.446 0.626 0.382 0.013 


642 -0.141 0.036 -0.131 0.837 -0.078 -0.097 


356 0.766 -0.055 0.020 0.127 0.907 -0.110 


127 0.708 -0.013 -0.331 -0.007 0.507 0.342 


Temperature difference average 
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差 均 值 作 自 


年 份 Year 
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0.001 -0.022 0.006 -0.002 0.954 


系数 Coefficient 


3 1992 一 2010 年 江西 省 温 湿 系数 滑动 均值 森林 病虫害 发 生 趋势 与 模型 拟 合 
Fig.3 Model fitting and trend of hydrothermal coefficient sliding average (HCSA), forest diseases and pest of Jiangxi during 1992 to 2010 


从 以 上 主 成 分 中 ,选择 荷载 最 大 的 温 湿 系数 请 动 均值 .相对 湿度 滑动 均值 .日照 时 数 变 率 、 降 水 量变 率 、 温 


变量 ,森林 病 虫 灾害 发 生 作 因 变 量 继续 进行 逐步 回归 ,获得 3 个 最 优 回归 方程 
Yrsg =3.582 x 105-7.750 x 105 X(R2=0.779，P<0.001 ) 


了 十 二 一 0.375X 10°X+2.95x10°( R* =0.791,F=49.335 ,P<0.001) 


Ya =—1.375x10°X+6.321x10°( R* =0.567,F=16.996,P=0.001) 

相似 地 ,以 上 X 均 为 温 湿 系 数 滑动 均值 ,三 个 方程 的 线性 拟 合 度 分 别 为 77.9% 、79.1% 与 56.7% ,模型 均 具 
预测 能 力 ,而 将 温 湿 系数 滑动 均值 回 代 后 ,得 知 其 平均 预测 准确 率 为 66.2% .68.6% .47.9% (图 3)。 因 而 ,与 
其 他 气候 要 素 相 比 , 温 湿 系 数 滑动 均值 可 能 对 区 内 和 森林 病 虫 灾害 发 生 及 趋势 形成 起 着 最 关键 的 影 


在 研究 区 内 呈 ! 


T 


3.8 


响 , 但 二 者 
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; 著 负 向 关系 ,前 者 在 过 去 几 十 年 间 显著 下 降 , 而 森林 病 虫 灾害 发 生 面积 则 显著 增加 (图 3)。 
2.2 ”江西 省 1961 一 2010 年 温 湿 系数 时 空 变化 特征 
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2.2.1 温 湿 系数 变化 倾向 率 与 周期 性 

研究 区 1961 一 2010 年 间 的 温 湿 系 数 滑动 均值 为 4.517 ,总 体 呈 显著 下 降 趋 势 ( 尼 = 0.5827，P<0.001 ,图 
4A) ,每 10 年 的 线性 倾向 率 为 -0.075, 且 正 、 负 距 平 转折 点 在 1993 年 (图 4A) 。 小 波 方差 图 (图 4B) 显 示 : 过 
去 50 年 间 区 内 温 湿 系数 变化 有 2 个 周期 ,小 周期 为 2a, 大 周期 为 29a, 而 以 后 者 为 主 。 小 波 实 部 变化 图 (图 
4C) 也 显示 29a 周期 振荡 最 显著 ,小 波 能 量 图 (图 4D) 最 大 值 1.4 也 出 现在 29a 这 一 时 间 尺 度 位 置 上 。 
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4 江西 省 1961 一 2010 年 温 湿 系数 滑动 均值 时 间 序列 变化 趋势 (A, * * P<0.001) 与 周期 性 特征 (B 小 波 方 差 图 \C 小 波 实 部 变化 图 .D 小 
波 能 量 图 ) 
Fig.4 ”Variation ，tendency and periodicity of hydrothermal coefficient sliding means during 1961 to 2010 in Jiangxi Province 


2.2.2 温 湿 系 数 变化 倾向 率 及 其 空间 分 布 特征 

江西 省 各 区 县 温 湿 系数 差异 较 大 ,区 域 性 变化 特征 显著 (图 5) 。 例 如 , 温 湿 系数 滑动 均值 总 体 上 沿 闽 西 
北 ( 最 高 5.59 ) 向 竟 东 南 (最 低 3.84) 递减 ,在 赣州 出 现 明 显 低 值 区 (图 5A)。 另 外 ,各 区 县 50a 温 湿 系数 滑动 
均值 线性 变化 倾向 率 或 变化 程度 则 沿 赣 西 .六 南 部 向 竟 东 . 议 北 部 递增 ,而 在 赣州 . 吉 泰 盆地 的 变化 程度 最 低 。 
由 此 说 明 ,过 去 50 年 间 ,闽南 高 温 或 干旱 更 重 , 而 竟 北 变动 性 更 大 。 


3 结论 与 讨论 


近 50 年 来 ,气候 变化 加 剧 .全球 趋 于 变 暖 "' 当 .陆地 表面 趋向 干旱 化 :5 ,各 类 森林 灾害 频繁 发 生 '“” ,损失 
加 重 '* 中 。 以 森林 生物 灾害 为 例 ,其 发 生态 势 近 几 十 年 来 在 我 国 呈 加 重 趋势 ,损失 也 随 之 增加 中 。 例 如 , 研 
究 区 江西 省 年 均 气 温 、 冬 季 均 温 50 年 来 (1961 一 2010) 分 别 以 0.16%7Z10al71 .0.27%Z10a 速率 升 高 ,已 分 别 升 
高 0.87°%C 、1.38%C ;然而 ,其 年 均 相 对 湿度 则 以 -0.45%/10a 的 速率 下 降 。 并 且 , 短 时 间 (1992 一 2010 年 ) 来 看 ， 
其 年 均 温 滑动 均值 夏季 均 温 滑动 均值 等 也 在 显著 升 高 ,而 年 降水 量 滑动 均值 相对 湿度 均值 相对 湿度 滑动 
均值 等 则 在 显著 下 降 , 这 些 足 以 说 明 研 究 区 气候 变 暖 及 和 干旱 化 显著 。 而 与 此 同时 ,区 内 森林 病 虫 灾害 发 生 也 
不 断 加 重 , 表 现在 灾害 面积 扩大 爆发 频繁 .一些 次 生性 或 偶发 性 害虫 局 部 爆发 增多 等 方面 。 据 范 霍 夫 定 律 ， 
气候 变 暖 驱 使 生物 发 育 速率 加 快 . 生 殖 力 增强 .存活 率 尤 其 是 越冬 存活 率 提 高 ,以 及 生物 发 生 始 见 期 . 迁 飞 期 、 
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N 

入 
3.8422 一 4.0431 -0.1700 一 -0.1523 
4.0431 一 4.1603 ma -0.1523—-0.1327 
Bl 4.1603 一 4.2286 上 -0.1327—-0.1149 
国 4.2286—4.2684 ma -0.1149—-0.0989 
4.2684—4.3367 a -0.0989—-0.0843 
ml 4.3367—4.4538 a -0.0843 一 一 0.0712 
有 4.4538 一 4.6548 a -0.0712 一 -0.05393 
ma 4.6548—4.9995 -0.0593 一 -0.0486 
ma 4.9995—4.5908 -0.0486 一 -0.0388 
ml 5.5908—4.6052 -0.0388 一 -0.0300 


图 5 江西 省 各 区 县 1961 一 2010 年 间 温 湿 系数 滑动 均值 (A) 及 其 变化 倾向 率 (B) 分 布 特征 
Fig.5 Spatial variation of hydrothermal coefficient sliding average (A) and change trend rate (B) 


高 峰 期 提前 ,年 生育 期 延长 .发 生 代数 增加 等 "1 ,从 而 促进 害虫 发 生 。 而 气候 变 暖 也 能 驱使 生物 适 生 区 加 宽 ， 
为 次 生性 或 偶发 性 害虫 发 生 提供 有 利 机 会 。 本 次 研究 支持 气候 变化 、 变 暖 驱使 森林 病 虫 灾 害 加 重 的 观点 ,如 
歼 淑 平等 认为 大 范围 流行 性 毁灭 性 病虫害 的 发 生发 展 \ 流 行 都 几乎 与 气候 变化 有 关 ,Dukes 等 '” 则 强调 
温度 上 升 3 一 5%C 可 加 重 球 蚜 .天幕 毛虫 蜜 环 菌 根 腐 病 和 山 毛 样 烂 皮 病 对 欧 鼠 李 、 南 蛇 苹 等 植物 的 侵害 ,从 而 
破坏 森林 结构 成分。 事实 上 ,本 次 所 选 24 个 气候 要 素 有 16 (与 病害 发 生 ) 或 17( 与 病虫害 、 虫 害 发 生 ) 个 显 
车 相关 要 素 ,夏季 均 温 滑动 均值 等 为 最 大 正 相 关 , 温 湿 系 数 滑动 均值 为 最 大 负 相 关 要 素 ,表明 过 去 儿 十 年 的 温 
度 上 升 (尤其 是 夏季 高 温 ) 或 环境 干旱 促使 着 森林 病 虫 灾害 发 生 加 重 ,这 与 一 些 研 究 认为 年 均 气 温 升 高 是 导 
致 森林 病 虫 灾害 加 重 的 重要 因素 观点 一 致 。 但 是 ,部 分 研究 倾向 于 支持 暖冬 的 重要 驱动 作用 5 ,如 1 月 份 暖 
温 驱使 杨 干 象 .美国 白 蛾 . 松 纵 坑 切 梢 小 索 等 森林 害虫 灾害 加 重 5 2 ; 而 另 一 部 分 则 倾向 于 支持 夏季 高 温 的 重 
要 驱动 作用 ,如 6 月 高 温 是 导致 洪都拉斯 南方 松 小 圳 虫 ( Dendroclomnas frontalis ) 爆发 的 重要 因素 之 一 5 ,本 研 
究 也 支持 森林 病 虫 灾 害 发 生 更 易 受 夏季 高 温 驱 动 。 

然而 ,气候 变 暧 与 环境 干旱 的 双重 或 联合 作用 可 能 较 其 单一 作用 更 为 强烈 。 在 分 析 各 气候 要 素 对 研究 区 
和 猴 林 病 虫 灾害 发 生 的 影响 中 , 先 将 所 有 气候 要 素 与 森林 病 虫 灾害 发 生 面积 进行 相关 分 析 , 往 出 与 之 显著 相关 
者 (16 或 17 个 要 素 ) 进 行 主 成 分 分 析 , 得 到 4 或 5 个 主 成 分 ,从 而 降低 维度 ,减少 共 线 性 干扰 。 在 获得 的 主 成 
分 中 ,发现 代表 温度 或 温 湿 度 联合 变量 特征 的 第 一 主 成 分 贡献 率 最 大 (病虫害 与 虫害 :41.43% ,病害 :42% ) ， 
该 结果 支持 温度 对 和 森林 病 虫 灾害 发 生 具 重要 影响 的 观点 ,而 又 突显 出 温 湿度 联合 作用 的 重要 作用 。 事 实 上 ， 
气候 变化 对 病 虫 灾害 发 生 影响 复杂 ,因为 单个 气候 要 素 的 变化 常会 引起 其 他 要 素 相 应 地 变化 ,如 生物 的 适宜 
温度 范围 可 能 因 湿 度 变 化 而 偏 移 , 反 之 亦 然 , 所 以 即便 是 不 同 的 温 湿 度 组 合 也 可 能 产生 出 近似 的 生物 发 生效 
应 55。 例如 ,落叶 松 钙 叶 蜂 的 发 生 就 同时 受到 温度 .湿度 .降水 3 个 气候 要 素 的 联合 作用 '”" ,同样 地 ,湘西 南 
马尾 松 毛虫 的 发 生 也 不 是 单个 气候 因子 所 左右 ,因为 山区 温度 .湿度 和 降水 都 能 对 其 发 生产 生 影响 。 尽 
管 ,从 前 的 一 些 研究 及 本 次 研究 结果 均 支 持 环 境 温度 变化 是 影响 森林 病 虫 灾害 的 主要 因子 的 观点 。 但 是 ,不 
容 忽视 的 是 ,环境 湿度 的 变化 也 可 能 参与 其 中 或 以 联合 效应 发 挥 作用 。 事 实 上 ,环境 干旱 对 一 些 害虫 发 生 直 
接 有 利 , 如 小 陇 山林 区 的 葡萄 长 须 卷 峨 发 生 与 危害 与 当地 温暖 干燥 的 气候 环境 密切 相关 "| ,意大利 阿尔 卑 斯 
山地 区 的 云 杉 叶 蜂 易 因 干旱 胁迫 或 在 气候 干 热 条 件 下 爆发 成 灾 '* ,而 锈 色 粒 肩 天 牛 也 多 在 降水 少 、 温 湿 系 数 
低 的 年 份 爆发 成 灾 '”]。 
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实际 上 ,在 预测 模型 构建 中 ,也 发 现 温 湿 系数 请 动 均 值 这 一 综合 要 素 表现 最 好 。 例 如 ,经 主 成 分 分 析 后 ， 
选择 各 主 成 分 中 载荷 值 最 大 者 为 自 变量 、 和 森林 病 虫 灾害 发 生 面积 作 因 变量 ,以 逐步 回归 法 构建 获得 的 最 优 回 
归 预 测 模型 , 均 是 基于 温 湿 系 数 滑 动 均值 (X) 所 构建 (Ys = 3.582 x 10 -7.750 x 10 X;Yas =-6.375x10 X+ 
2.95x10 ;Yi = 一 1.375x10 外 +6.321x10 ), 且 与 用 主 成 分 回归 法 所 获 模 型 (如 Yrs = -3377.397Y,+ 
1352975.004) 相 比 ,前 法 所 构建 的 模型 拟 合 度 更 高 (如 拟 合 度 Yas = -6.375x10*X+2.95x10 为 79.1% ,而 
Yung = 一 3377.397Y,+1352975.004 为 62.7% ) 。 作 为 年 均 相 对 湿度 与 年 均 温度 的 比值 , 温 湿 系 数 是 反应 温 湿度 
联合 关系 与 作用 的 综合 性 气候 要 素 , 年 均 相 对 湿度 下 降 或 年 均 温度 上 升 ,或 者 两 种 变化 同时 发 生 ,导致 温 湿 系 
数 下 降 。 研 究 表明 ,江西 省 过 去 50 年 年 均 温 滑动 均值 明显 上 升 ,而 相对 湿度 均值 与 相对 湿度 滑动 均值 则 明显 
下 降 , 温 湿 系数 滑动 均值 因此 呈 显 著 下 降 趋 势 ,为 一 些 森 林 病 虫 灾害 发 生 创 造 了 有 利 条 件 交 2 3 ,从 而 驱使 
着 森林 病 虫 灾害 在 过 去 几 十 年 间 不 断 加 重 。 尤 其 , 随 着 全 球 在 20 世纪 90 年 代 进 入 近 百 年 来 最 暖 的 10al?3) ， 
温 湿 系数 下 降 也 最 为 明显 ,如 其 滑动 均值 即 在 1993 年 转 为 负 距 平 ,驱使 着 森林 病 虫 灾害 深度 加 重 。 并 且 , 研 
究 区 温 湿 系 数 滑动 均值 具 29a 大 周期 特征 , 故 若 以 1993 年 为 切入 点 , 则 可 预 判 区 内 森林 病 虫 灾害 加 重 趋势 至 
少 维持 到 2022 年 。 男 外 ,研究 区 温 湿 系数 滑动 均值 具 区 域 分 布 特征 ,如 温 湿 系 数 滑动 均值 总 体 沿 座 西北 向 六 
东南 递减 ,并 在 闵 州 市 及 其 周边 形成 一 个 低 值 区 (图 5A)。 尽 管 座 州 . 吉 泰 倪 地 等 地 为 江西 省 森林 病 虫 重 灾 
区 ,但 据 温 湿 系 数 滑动 均值 线性 变化 倾向 率 呈 沿 着 诸 西 、. 闭 南 部 向 赣 东 、. 议 北部 递增 ,可 预 判 随 着 气候 变 暖 或 
环境 干旱 , 竟 东 北 森 林 因 病 虫 灾害 而 遭受 破坏 的 加 重 进度 将 较 其 他 区 域 更 快 ,换言之 , 轻 灾 区 的 未 来 森林 病 虫 
灾害 变动 性 可 能 较 重 灾区 更 大 。 


致谢 :本 文 地 面 气候 资料 日 值 数 据 来 自 中 国 气 象 数据 中 心 ,江西 农业 大 学 林学 院 研究 生 曾 城 做 了 大 量 的 数据 
录入 工作 ,特此 致谢 。 
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